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Deriving Mathisson - Papapetrou equations from relativistic 

pseudomechanics 

R.R. Lompay* 


Abstract 

It is shown that the equations of motion of a test point particle with spin in a given gravitational held, so called 
Mathisson - Papapetrou equations, can be derived from Euler - Lagrange equations of the relativistic pseudome¬ 
chanics - relativistic mechanics, which side by side uses the conventional (commuting) and Grassmannian (anti¬ 
commuting) variables. In this approach the known difficulties of the Mathisson - Papapetrou equations, namely, 
the problem of the choice of supplementary conditions and the problem of higher derivatives are not appear. 

PACS numbers: 04.20.-q, 12.60.Jv 

In ref. |T| dealing with a description of classical relativistic point particle with spin was suggested a model based 
on action functional 

e 2 

W[z, ip; 9 1 , 2 } = 1 J @9 

81 

Here c is light velocity in vacuum; Ti - Planck constant; m - bare mass of the particle; n - dimensionless coupling 
constant; 9 - an arbitrary parameter which numerates the points of the world line of the particle and has length 
dimension; z = = (z°, z 1 , z 2 , z 3 ) are pseudocartesian coordinates of the particle in Minkowski space; = {7 ^ A b } 

Dirac matrices which satisfy anticommutation relations 

T Y + tY = -2 rf v I, (2) 

where 

{Vuv} = {rf v } = diag(-l, 1,1,1) (3) 

is Minkowski space metrics; I = {& A b} ~ unit 4 x 4-matrix; ip = {ip A } and ip = {ip B } = ip ^7° are Dirac spinor and 
conjugated to it spinor, accordingly, which are defined on the world line of the paricle and considered as Grassmannian 
variables: 

iP a iP b = -ip B ip A , ip A ip B = - ip B ip A , ip A ip B = -ip B ip A ■ (4) 

@9 = [— r] a pz' a (9)z' ,3 (9)] 1 / 2 d9 and D = [— r\ a pz la (9)z l P(9')]~ x l 2 dld9 are, accordingly, reparametrization-invariant 
(RI) measure and Rl-derivative - the objects which allow to fulfill the calculations not depending on the special choice 
of the parameter 9 (see ref. 0); z' a (0) = dz a {0)/d,0. 

The model (1) correctly describes the kinematical connection between orbital and spin parts of the total 4-angular 
momentum J po = M pa + S pa and gives also true results for Poisson brackets for dynamical variables in corresponding 
Hamilton formalism. 

In this work we show that in the case of interaction such particle with an external gravitational field the equation 
of motion which follow from corresponding generalization of the action functional (1), leads to Mathisson Papapetrou 
equation (MPE) {30. 

We introduce an interaction of the particle with gravitational field with minimal coupling by covariantization of 
action functional (1). Metrics (3) and 7-matrices (2) of Minkowski space, which we shall denote now as ?7( a )(/3) and 
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—me 2 + {ip {9) Dip (9) — Dip{9)ip{9)) 


' 1 -mc 2 ip{9)^ a ip{9)Dz a {9) 


(1) 
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7 ' a ', we replace by metrics g pl ,(x) and 7 -matrices ^(x) of an arbitrary pseudorimannian space. And Rl-deriative D 
we replace by absolute Rl-deriative V. Then the action functional (1) take the form 

W[z,ip,ip,g-,6 1, 2 ] = 

@2 _ _ _ 

= i J 39 [-me 2 + i^(ip(9)Vip(9) - Vip(9)ip(9)) + %mc 2 ip(9)y a (z(9))ip(9)Dz a (9)} (5) 

V 

= w 3 + w 2 + w 3 , 

where now 39 = [—g a0 (z(9))z' a (9)z , d(9)] 1 / 2 d9; D = [—g a0 (z(9))z' a (9)z ,/3 (9)]~ 1 ^ 2 d/d9-, 7^(2;) = e^ a )^ (x) and 

e (a) M ( x ) form a fierbein: 

g (a){0) e {a f{x)e {0) v {x) = g^(x); g^ I/ (x)e ia f(x)e w ,/ (x) = ? 7 («)(/ 3 ), ( 6 ) 

so that 

7* 1 (x)y v (x) + y , '(x)y ft (x) = —2 g^ v (x)I. (7) 

Absolute Rl-derivative V acting on the spinors is defined as 

VV>(0) = D^(9) + ^y(z(9))iP(9)Dz v (x)-, 

Vip(9) = Di!>{9) - ip(9)VL v (z(9))Dz v (x), K ’ 

where fl„(:r) are spin connection coeffitients: 

^•v(x) = — ^ix( p ^ a ^ a ){x)e^ \(x)\ = e(p) fi (^)(VAe( tr )j 1 (i))e( 0 ) ( x ). (9) 

And the action on the vector is defined as 

VDz^(9) = D 2 z^(9) + T )1 Xv {z{9))Dz x {9)Dz v {9), (10) 

where T lx \ v (x) are Kristoffcl symbols: 

r M Ai/(*) = g^ a {x)]^(d\g av (x) + d v g a \(x) - d a g Xv (x)). ( 11 ) 

Let as find the variations of the action functional ( 8 ), caused by varying of quantities z(9), ip (9), if>(9): 

5W[z, ip, ip, g ; 6*1,2] = W[z + 5z,ip + 5ip,ip + Sip, g; 6*1,2] - W[z, ip, ip, g; 6*1,2] = 5 Z W + 5^W + SjW. (12) 

8 z W[z, ip, ip, g; 6 * 1)2 ] = W[z + Sz, ip, ip, g; 0 i, 2 ] - W[z, ip, ip, g\ 6*1,2] = 5 Z W\ + S Z W 2 + S Z W 3 ; ( 13 ) 


q 2 

5 Z W\ = —me f 5(39 ) = —me f d95 z ( x /-g a0 (z)z ,a z't 3 ) 


O 2 


= n rJd9 \-g p(J (z)z' p z'°\- x/2 (d p g af3 8z»z' a z'P + 2 g ap z' a (8z»y) 


= tf J 39 [ d fl g a0 Dz a Dz /S 5z» + 2 g ail Dz a D(5z IJ )\ 

0i 

= if S f 39 [d p g a0 Dz a Dzf } - 2 D(g Qp Dz a )\ 5z» + h °^ Y 


O 2 


= -me J 39 [g ap D 2 z a + y(d a g 0p + d 0 g afl - d lz g a p)Dz a Dz 0 ] 5z fl 

0i 

= J39 [-V(meg Xp Dz>)]5z»+ b ° Undary • 

0 i 


terms ’ 


boundary 

terms 


(14) 
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S z W 2 = *j J dO ipS z {ip'+ fl a {z)ipz la ) — S z {ip - ip£l a {z)z' a )ip =ih J d9ip5 z {fl a {z)z ,a ))ip 

0i L J 0i 

02 _ _ 

= ih J d9 \ipd l _ l fl a {z)ipz' a Sz IJ ' + ipQ fl {z)ip{Sz IJ 'y] 

0i 

@2 _ _ 

= ih J 99 \ipd t _ l fl a {z)ipDz a Sz IJ ’ + ip^l /J ,{z)ipD{5z IJ ’)] 

0i 

= m j m $(dM^Dz a - Dm^)] sz» + bo t ^ ry 

@2 _ _ _ 

= ih f 99 ^p{d^,Cl a — d a fl^,)ipDz a — Dipfl^ip — ipfl^D'ip] Sz M 

0i 

@2 _ 

= ih f 99 \$(diSl a — 3 a f2 M + — £l a Qp)ipDz a — 

0i 


boundary 

terms 


{Dip - tpSl a Dz a )Sl lt ip - ipn^Dip + ipQ. a Dz a )\ Sz* + boundary 

termb 

Jm [±(h^ S P°iP)R p ^ a Dz a - ih(V^%iP+ ^X7iP)\ 6z» + bo t ^ ry ; 


S Z W 3 = %mc J d9ip5 z (j a (z)z' a ))ip = pnc f 99 ^{d^ajipDz 01 - D(ipj^ip)] 5z» ■' boundaiy 

0i 0i 

02 


terms 


f he [?(^7« - I* a „7 fiWDz* - {D(ip 7^) - rV(</>7 vip)Dz a )] 5z» + b °^^ 


= / 99 [^mc’tp{'ja^ l i - £L^ a )ipDz a - V(§racV’7 i uV’)] $ zM + 

0i 

In obtaining the last equalities in (15) and (16) we used the identities 


boundary 

terms 


7" HyiHa Qa^fi — 2! SP 


s p ° = -{Yi a -1”Y)- 


^pHot — d^c ^^a/i'Y/3 “b Qfj.'Ja 7 a^fx — 0* 

Furthermore, 

5$W[z, ip, ip, g; 6b, 2 ] = W[z, ip + Sip, ip, g; 0 1)2 ] - W[z, ip, ip, g; 0i )2 ] = S^W 2 + S^W 3 -, 

O 2 &2 

S^W 2 = f / 99 [PpS^Dip + n a ipDz a ) - V^V>] =ff 99 $D{6ip) +i^ a SipDz a - V^Sip] 

0i 0i 

= § / ^ [-(^ - W a Dz“) - V#V>] <¥ + b ° t e^ a s iy 


= f 99 —ih\7ip Sip ■ 


boundary 

terms 


S^W 3 = / 99 [P;mcip'y a Dz a ] Sip. 

0i 

Analogously 

^W[z, i>, 0; 01 ,2] = W|>, V>, V> + 0; 01,2] - W[z, g; 01,2] = <%W 2 + < 5 ^W 3 ; 


02 


S^W 2 = / 99 Sip [ih\7ip\ 


boundary 
terms ’ 


02 _ 

fcIF 3 = / 00<ty> [fTOC7 a iAD 2 a ] . 

0i 


( 15 ) 


(16) 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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Taking, as usually, that on the boundary of integration 9 = 9 \ 2 , the variations vanish, 8 z \ t 2 
using the formulae (12) (16), (19) (24), we obtain the following equations of the motion 


AW 

A z p 


-V 


mcg X p(Dz x + ^ip 7V) + ^( hips pa ip)R P ap V Dz v - 


— (ihVip — '^mc'ip'y a Dz a ) VL^ip — ipQ.p[iKShp + ^ mc"f a 4’Dz a ) = 0; 


A R W 
Aip 

A L W 

At/) 


— (ihVip — ^mc^uD z a ) = 0; 
ihVip + —mc"fa4’Dz°‘ = 0. 


Substituting the eqs. (26), (27) into eq. (25), we reduce the last system to the form 


Slpl,2 = Sip i2 = 0, 

(25) 


(26) 

(27) 


f V [mc{Dz p + § ip^ip)\ = ±R» V p <T Dz v {hips pri xp)-, 

) V^=i K (^) 7 a ^Dz a -, 

{ ViP = -^K(^)^ a Dz a . 

Using the notation 

pp =m [Dz P + ^lp^{z)lp^ , Spa = lpSpa{z)lp 


(28) 


(29) 


for the expressions which are the general covariant generalization of the expressions for 4-momentum and 4-spin of 
the particles in pseudoeuclidian space-time PQ , we can write the first eq. in (28) also in the form 

VP 11 = ^Dz^. (30) 

With the help of eqs. (28),(29) it is easy to obtain the equations 

VSpa = -(DZpPa - DZaPp ), (31) 

As consequences of these equation one obtain the equations 

VSpa + D Z pDz"S„a + Dz a Dz v S pv = 0, (32) 

and 


Dz v X/S„p = (<5% + Dz v Dzp)P, (33) 
Using the eqs. (33), one can represent the 4-momentum in the form 

PP = —Dz p Dz u P 1 ' + (5 V p + Dz u D Z p)P l/ = MDz p + Dz^VS^, M =-Dz v P v (34) 

Substituting the decomposition (34) into the eq. (30), one obtain 

V(MDz p + Dz u VS„ p ) = ^RP vp aDz l 'S pa . (35) 

The system of eqs. (32), (35) 

f VS pa +Dz p Dz v S va + Dz )T Dz"S pv = 0; ( . 

\ V(MDz p + Dz^S^) = ±RP vp aDz v S pG 1 j 

at first were obtain in papers El El and now is known as Mathisson Papapetrou equations. The first equation of the 
system is called the equation of motion of the spin and second one - the equation of motion of the particle in space. 

In original papers El El these equations were derived halfphenomenologically by integrating of multiple 3-momentum 
of a matter distributing in space with posterior reducing it to a point and covariantization the obtaining results. In 
our consideration the MPE are obtained as exact consequences of the variation Euler - Lagrange equations for abhoc 
point generally covariant physical system. This approach is automatically free from the well known difficulties MPE, 
namely. 
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The first is the problem of choice of supplementary conditions. If the eqs. (36) are considered as the closed system 
for determining z M (0) and S pa {8 ), then it is easy to see, that the number of equations is less then the number of 
unknown functions. The first eq. (36) determines not the total \7S pa , but only its space projection 

KKVS pu = (,5 \ + Dz»Dz p )(5\ + Dz^Dz a )VS pu = VS pCT + Dz p Dz v S viy + Dz a Dz v S pv . (37) 

So, the system MPE is not fully determined. For full determination one puts on S pa some supplementary conditions, 
such as 

Corinaldesi-Papapetrou-Pirani condition 13 El S^Dz v = 0; (38) 

Tulczyjew condition [2] P v = 0 (39) 

and some others. The problem is that there is not essential reason to prefer one or other supplementary condition. 
In our approach based on action functional (5) this problem doesn’t arises. Indeed, the eqs. (36) are only a sequence 
of the basic system of eqs. (28) which fully determines all fundamental variables z(8), ip(9), if)(9) when initial data 
zq = z(9o),ipo = i/j(9o),ip 0 = are given. By known fundamental variables all other dynamical variables are 

determined unique. So, for example, the 4-spin is given by the second formula in (29). 

The second is the problem of higher derivatives. As it is easy to see, the left-hand side of the second eq. in (36) 
involves the second derivative of S pa , which is not determined, since, as were noted, the first eq. in (36) fixes only 
the space part of S/Spa- When the supplementary condition (38) is imposed, the expression in parentheses in the 
second eq. in (36) can be rewritten as ( MDz^ — V(Dz' y )S u tJ '). Then the left-hand side of the second eq. in (36) will 
include the third derivative of coordinates of the particle, but this contradicts to traditional form of mechanics. In our 
approach this problem is not arising, too. In accordance with the equations given above, the expression in parentheses 
in second eq. in (36) is equal to m(Dz M and higher derivatives are absent. 
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